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Abstract 

We study under-barrier tunneling for a pair of energetically bound bosonic atoms in an optical 
lattice with a barrier. We identify conditions under which this exotic molecule tunnels as a point 
particle with the coordinate given by the bound pair center of mass and discuss the atomic co- 
tunneling beyond this regime. In particular, we quantitatively analyze resonantly enhanced co- 
tunneling, where two interacting atoms penetrate the barrier with higher probability than a single 
atom. 
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I. INTRODUCTION 



The phenomenon of under-barrier tunnehng is one of the most exciting predictions of 
quantum mechanics, which does not fit the classical picture of the world. As a bright 
example of a pure quantum effect, it is considered in all text books on quantum mechanics, 
which might make an impression that under-barrier tunneling has been exhaustively studied 
already in early days of quantum mechanics. However, the theoretical description of the 
under-barrier tunneling is simple only for a point particle in one dimension. If we have 

a composite object, the problem of tunneling becomes very involved. Systematically this 

■ i 

problem was first addressed in nuclear physics (see [li] and references therein). It was found 
that the tunneling probability for the composite object may considerably differ from that 
for the point particle of the same mass. 

In the present work we revisit the problem of the under-barrier tunneling for a composite 
object which has become experimentally available only recently - the pair of interacting 
bosonic atoms in an optical lattice, where two atoms stay close to each other due to the 
energy constrain ^j. Note that such atomic pairs exist for both attractive and repulsive 
inter-atomic interactions if one satisfies the necessary condition that the interaction energy 
\U\ exceeds the single-atom tunneling energy J. Following Ref. we shall refer to this 
exotic two-atom molecule as the bound pair. It admits a simple yet rigorous description ^ 
^, which greatly facilitates theoretical studies of different aspects of the composite object 
tunneling. 

Our other motivation for studying the under-barrier tunneling of the boundpair are 
problems of macroscopic tunneling of a Bose-Einstein condensate of cold atoms 7H1C|. In 
particular, the recent paper [10] analyzes numerically the under-barrier tunneling of a bright 
soliton consisting of ^ 1 condensed atoms. Since the bound pair can be considered as 



the bright soliton with N = 2 atoms 



a 



rigorous analysis of the bound-pair tunneling 



may shed additional light on the phenomenon of macroscopic tunneling. 

The structure of the paper is as follows. In Sec. [Tllwe recall the main results on eigenstates 
of the bound pair in the absence of an external potential and introduce a simple two-state 
model, which suffices to describe the mobility of the bound pair. Sec. Hill is devoted to 
the under-barrier tunneling within the framework of the two-state model. We identify the 
conditions under which the bound pair tunnels as the point particle with the coordinate given 
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by the center of mass of the bound pair and uncover the effect of resonant tunnehng, which 
is entirely due to internal degrees of freedom of a composite object. The main drawback 
of the two-state model is that it neglects the dissociation process where the barrier breaks 
the pair into two unbound atoms. For this reason, in Sec. IIVI we simulate the tunneling 
process numerically without using any approximations. We summarize our findings in the 
concluding Sec. |Vl 



II. TWO-STATE AND ONE-STATE MODELS 



As mentioned in the introductory section, the strongly interacting bosons in a lattice 
form bound pairs, where two bosons occupy the same site. Such a pair can move across the 
lattice by virtually breaking the bond 12| and, thus, is a composite object with well defined 
kinetic energy. The dispersion relation E{k,) for the bound pair can be easily calculated 
numerically by diagonalizing the Hamiltonian of the Bose-Hubbard model [see Eq. ([T]) below] 
with N = 2 particles. For the purpose of future references Fig. [T](a) shows the result of this 
diagonalization for a lattice comprising 11 sites, where we additionally parametrize the Bose- 
Hubbard Hamiltonian by the Peierls phase 9: aj+ifl; — exp(z^). In Fig. [U^a) the bound 
pair of two bosons is associated with the lower band, while the upper band is the spectrum 
of two hard-core bosons. The problem can be also solved analytically, either exactly or by 
using a perturbative approach. In the rest of this section we discuss the dispersion relation 
E{k,) and the eigenstates of the bound pair within the perturbative approach, which better 
fits our aims of studying the tunneling process. 

To facilitate the theoretical analysis it is convenient to consider the Bose-Hubbard model 
which also includes interactions in neighboring sites: 

Hbh = -^^{al+iai + h.c)j + ^^ni{ni - 1) + Ui^hi+ihi . (1) 

Inclusion of the latter term explicitly introduces excited states of the bound pair, see 
ig. [H^b). We note that we do not assign the interaction constant Ui any physical meaning 



131] and let eventually Ui tend to zero. The perturbative approach to the energy spectrum 



of the bound pair essentially amounts truncation the Hilbert space of the operator ([T]) to 
the subspace which includes only the Fock states where two bosons occupy either the same 
site or two neighboring sites. Then, denoting by \E'P^ the probability amplitude to find two 
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bosons at the site I and by to find them at sites / and / + 1, the eigenvalue equation for 
the bound pair takes the form (ifo^)z = -E'^h where 



V2J 



1 1 

10/ I 1 / loo 




(2) 



The solutions of this eigenvalue problem are plane waves, \E'; = C(K)e*'^', where the vector 
C{k) satisfies the following 2x2 eigenvalue equation: 



A 



-J{l + e"')/V2 



V 



-J(l + e 
From ([3]) we have 

E{k) = Uo 
In what follows we refer to Eqs 







C = EC, A = \Uo-Ui\ 



(3) 



A 



± 



A 



2J2 



cos^ 



(4) 



2 VV2/ V2, 
as the two-state model. 
It is worth stressing that the two-state model provides only an approximation to the 
exact eigenstates of the bound pair. How good this approximation is depends on the system 
parameters. In general, the two-state model is a good approximation for both the ground 
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FIG. 1: The band spectrum of the system ([T]) for J = 1 and {Uo,Ui) = —(2,0) (left panel) and 
{Uq,Ui) = —(5,3) (right panel). The lattice comprises 11 sites with periodic boundary conditions. 
The red lines are Eq. ([5]) and dH), respectively. 
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and excited bands if |f/i|, |?7o| ^ J and A ^ |f/o|. If is decreased, yet |?7o| ^ J, 
it is still a reasonable approximation for the ground band [minus sign in Eq. (j4])]. This 
also includes the case Ui = where the excited band is 'dissolved' into the spectrum of 
unbound bosons. In this case, according to our numerical analysis, \Uq\ should be at least 
twice as large as J. Then the admixture of the truncated Fock states (i.e., those belonging 
to the truncated subspace of the Hilbert space) to the exact ground state of the bound pair 
does not exceed 5 percents. For smaller \Uo\ the contribution of these Fock states cannot 
be neglected and one has to diagonalize the Hamiltonian ([1]) in the whole Hilbert space {4], 
which results in the dispersion relation 

E{k) = -^U^ + A.P cos2(/t/2) . (5) 

In Fig. [1] we plot the analytical results ([5]) and (|4]) by the red hues. 

If the band gap A ^ J the problem can be simplified further, resulting in the one-state 
model. The procedure goes as follows. First we restrict ourselves to the ground band and 
introduce the Wannier states of the bound pair by integrating its translationally invariant 
eigenstate over the quasimomentum in the first Brillouin zone: 

= r ^,e-*"'dK = r C(-)(/s:)e'(''-')MK . (6) 

J —TV J —7T 

The Wannier states (jS]) are localized functions with the center of gravity at the site /. (For 
example, for J = 1, [/g = -4 and Ui = we have |^o) ~ 0.157| . . . , 0, 1, 1, 0, 0, . . .) + 
0.975| . . . , 0, 0, 2, 0, 0, . . .) + 0.1571 . . . , 0, 0, 1, 1, . . . , ).) Next we calculate matrix elements 
of the Hamiltonian ([2]) for the Wannier states seperated by m sites: = (^z+ml-f^ol"^;)- 
This way we obtain the effective Hamiltonian where the bound pair is considered as a point 
particle: 

Heff = -lY.IrnE{blJl + h.C.) • (7) 

In the limit A ^ J we have Ji = 2J^/A and one can safely neglect the next to neighboring 
hopping. Obviously, this situation corresponds to the case where the dispersion relations 
(HllSll are approximated by the cosine function (which in practice requires A > 4J). 

III. UNDER-BARRIER TUNNELING FOR THE TWO-STATE MODEL 

The one-state model ([7j) introduced in the previous section gives us the reference frame 
in studying the under-barrier tunneling of the bound pair. To be certain we shall consider a 
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Gaussian barrier, e; = V exp{—P /2(t'^) , and a plane wave coming from minus infinity. Since 
the Gaussian barrier is well localized within the finite interval |/| < L ~ o" we can find the 
tunneling probability by using, for example, the transfer matrix method (see Appendix). 
Alternatively, one finds the tunneling probability by simulating the scattering process for a 
localized wave packet on the basis of the time-dependent Schrodinder equation: 



As the initial conditions for (|8]) it is convenient to choose a wide Gaussian with the given 
group velocity, ipi{t = 0) = G{1 — Iq) exp{iK,l). If the width of this initial packet is large 
enough, the result of time-dependent simulations practically coincides with that obtained 
on the basis of the stationary Schrodinder equation. 

A remark concerning the sign of the parameter V is in order. In what follows we consider 
both positive and negative V, i.e., potential barriers and wells. In fact, for a particle in a 
lattice the notions of 'barrier' and 'well' are equivalent to each other to some extent. This 
becomes especially clear in the case of neighboring hopping (1^ = /i5m,i), where scattering 
of the plane wave with the quasimomentum n on the well {V < 0) is equivalent to scattering 
of the plane wave with the quasimomentum k' = vr — k on the barrier {V > 0). [Note that 
for K = 7r/2 this implies a symmetric function Pt{V) for the tunneling probability.] Also, by 
considering both positive and negative V we cover the case of repulsive interactions as well, 
with the obvious substitution V ^ —V when the sign of interaction constants is changed. 

All said above about the one-state model is equally applied to the two-state model. Here, 
instead of (|8]), one deals with the Schrodinder equation 



where is a two-component vector and the Hamiltonian Hq is defined in Eq. ([2]). Our 
particular interest is the scattering of a plane wave transmitting in the ground energy band. 
As an example. Fig. [2] shows the tunneling probability for the plane wave with k, = tt/2 for 
3 different values of the parameter A. In the panes (b,c) of this figure we also plotted the 
tunneling probability obtained on the basis of the one-state model. 

A remarkable prediction of the two-state model as compared to the one-state model is 
the appearance of narrow transparency windows for negative V. Usually such windows are 



idtipi = -- 



(8) 



m 
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FIG. 2: Two-state model: Tunneling probability as the function of V for fixed k = tt/2 and 
A = 1,2,4. The scattering potential is e/ = 1/ exp(— P/2o"^), a = 0.65. The open circles are 
predictions of the one-state model. 

associated with resonant tunneling in multi-barrier structures. In our case (single barrier or 
well) we meet a different type of resonant tunneling, where the bound pair tunnels through 
the upper band or, more precisely, through a localized state of the bound pair in the excited 
state, as it is pictorially shown in Fig. |3l This interpretation of the enhanced tunneling is 
strongly supported by numerical simulations of the wave-packet dynamics, where we observe 
a temporal population of the upper band when the packet passes through the well. This is 
shown by the dashed line in Fig. HI where the upper curve is the total occupation probability 
of Fock states with two bosons in the same site and the middle curve is the total occupation 
probability of Fock states with two bosons in the neighboring sites. 

Number and widths of resonances seen in Pt{V) crucially depend on the system param- 
eters, in particular, on the width a of the potential well. If a is increased, we observe more 
resonances and they are narrower. The decrease of A makes resonances wider. It should be 
also mentioned that, for the currently considered Gaussian potential, the resonances appear 
in pairs, as seen in Fig. EJ^c), and the pair can merge into the single wide resonance, as it is 
the case depicted in Fig. |2]|^a-b). 
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FIG. 3: Pictorial presentation of the resonant tunneling. The letters (d-c) labels the ground and 
excited energy bands of the bound pair, (b) is the energy band of two bosons with one of them 
captured in the well, and (a) is for unbounded bosons. 

IV. THE FULL SYSTEM 

Next we discuss the degree of validity of the two-state model. Indeed, the two-state model 
neglects the coupling to the truncated Fock states, which are associated with unbound 
bosons. If this coupling is strong (as in the case t/i ~ 0, where the upper band of the 
two-state model is embedded into the energy band of hard-core bosons) it may essentially 
affect the tunneling process and even open new scattering channels where the bound pair 
dissociates. For this reason we simulate the tunneling process on the basis of the Bose- 
Hubbard model ([1]), i.e., without using any approximations. In this numerical experiment 
we propagate the wide Gaussian packet constructed from eigenstates of the bound pair for 
a time approximately twice longer than required for the packet to hit the potential barrier. 
Fig. [5] shows a typical result for [/q = —2, Ui = and V = —2. The figure depicts 
probabilities to find the system in the Fock state with one boson at the site / and the other 
one at the site m. (Note that m > / in the considered case of identical particles.) It is 
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FIG. 4: The total occupation probability of Fock states with two bosons in the same site (upper 
curves), Fock states with two bosons in neighboring sites (middle curves), and Fock states with two 
bosons separated by empty sites (lower curves). The dashed, dash-doted, and solid lines distinguish 
between the two-state model, the three-state model and the original system. Parameters are J = 1, 
(Uo, Ui) = -(4, 2) (hence A = 2),a = 0.65, and V = -2.2. 

seen in Fig. [5] that the initial packet splits into four packets where two of them, which are 
located at the main diagonal, are associated with the bound pair and the other two are the 
dissociated pair with one boson staying in the potential well. (For the chosen parameter 
the dissociation is energetically allowed because the boson in the well accumulates almost 
the whole bound energy.) Summing up probabilities for this four packets, which are well 
separated in the Fock space, we find the tunneling, reflection, and dissociation probabilities. 

To systematically study the effect of the Hilbert space truncation (i.e., the effect of 
unbound bosons) we fixed the parameter A = \Uo — Ui\ and vary the interaction energy 
Ui, where the limiting case corresponds to Ui = 0. The parameter k, which defines the 
group velocity of the incoming wave packet, is fixed to k = 7i/2 and the parameter a, 
which defines the width of the potential barrier, is a = 0.65 (then the scattering potential 
essentially comprises 3 lattice sites). The results of our numerical simulations are depicted 
in Fig. [6] by symbols, which are connected by the dotted line to guide eye. Open circles 



9 



-60 
-40 




40 
60 

-60 -40 -20 20 40 60 
I 



FIG. 5: Result of numerical simulations for the wave-packet dynamics of the original system for 
J = 1, Uq = —2, Ui = 0, (7 = 0.65, V = —2. The figure encodes (by using the gray-scaled mapping) 
the probability to find two bosons at sites / and m. 




FIG. 6: Tunneling (open circles) and dissociation (asterisks) probabilities vs. V. The system 
parameters are k = it/2 J = 1, a = 0.65, and {Uo,Ui) = —(8,6) (top), —(4,2) (middle), and 
— (2,0) (bottom). The solid line is the prediction of the two-state model. 
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show the tunnehng probabihty and asterisks the dissociation probabihty. (The reflection 
probabihty, which is obviously given by = I — Pt — Pd, is not shown.) By inspection 
of the numerical data we can draw the following conclusions: (i) There are practically no 
deviations from the predictions of the one- and two-state models for the main transparency 
window around V = 0. Thus, here the bound pair tunnels just like a point particle; (ii) 
The system exhibits resonant tunneling at ~ A, as predicted by the two-state model. 
This is a clear manifestation of complex structure of the tunneling object; (iii) As compared 
to the two-state model the resonant tunneling is suppressed. It is interesting to note that 
the tunneling is suppressed independent of whether the dissociation channel is open or 
closed; (iv) If the dissociation channel is open, we observe strong back action of the resonant 
tunneling on the dissociation process, which manifests itself in a local deep in P^ = Pa{V). 

Let us discuss the suppression of the resonant tunneling in some more details. We found 
this suppression to be fairly reproduced if the Hilbert space of the two-state model is enlarged 
by including the Fock states, where two bosons are separated by one empty site. In the 
other words, instead of the two-state model one considers a three-state model. The dash- 
dotted lines in Fig. H] show population dynamics for the considered three families of Fock 
states. It is seen that it is practically coincides with that for the full system. Note that 
the third state remains practically unpopulated during the tunnehng process. Nevertheless, 
the presence of this third state appears to be important. This statement is also supported 
by the transfer matrix analysis, where the inclusion of this state considerably modifies the 
resonant tunneling (see Appendix). 

V. CONCLUSIONS 

We studied the tunneling of an interactively bound pair of two bosons in ID lattice 
through a narrow potential barrier/well. This system is, perhaps, the simplest composite 
object which can be created in the laboratory. We address the question under which condi- 
tions this composite object tunnels like a point particle. Loosely speaking, these conditions 
amount to the requirement that the microscopic interaction constant Uq entering the Bose- 
Hubbard Hamiltonian is larger than the hight of the potential barrier (the case of repulsive 
interactions) or the depth of the potential well (attractive interactions). If this condition is 
satisfied, the pair can be considered as a point particle. Note that the effective Hamiltonian 
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for this 'point particle' contains next to neighboring hopping, which is absent in the original 
Bose-Hubbard Hamiltonian. This takes into account the finite size of our composite object 
when it is treated point object 

HQ 

If the above condition is violated, one meets two phenomena which are entirely due to 
internal degrees of freedom of the composite object. These are (i) the resonant tunneling 
and (ii) the dissociation. Under conditions of resonant tunneling the bound pair can tunnel 
through the barrier, which the single boson cannot penetrate. (Here we refer to the case of 
repulsive interactions.) Assuming the shape of the barrier to be fixed, resonant tunneling 
takes place in a rather small region of the parameter space spanned by the quasimomentum k 
of the incoming plane wave (or group velocity of the incoming wave packet) and the height V 
of the potential barrier. On the contrary, the parameter region where the external potential 
breaks the pair by capturing one boson at the barrier is relatively large. For a generic form 
of the external potential and f/i = the region of resonant tunneling is usually embedded 
into the dissociation region and, thus, the resonant tunneling and dissociation coexist. We 
observed strong mutual infiuence of these processes which results in suppression of both the 
tunneling and dissociation. 

In the present work we calculated the tunnelin g an d dissociation probabilities by simulat- 
ing the wave-packet dynamics of the bound pair 15| where, as the physical object, we had 
in mind ultra-cold atoms in the ID optical lattice with the scattering potential created by 
additional laser beam crossing the lattice at the right angle. For this reason we considered 
a Gaussian shape of the potential barrier/well. It seems unlikely that one can obtain a 
compact analytical expression for the tunneling and dissociation amplitudes for this scat- 
tering potential. However, there are good prospects in developing the analytic theory for 
particular shapes of the barrier/well, which includes the impurity-like potential e/ = V5ifl 
and the box-like potential e/ = V{5i^q + 5; i). The advantage of the latter potential is that 
the well size and the size of the bound pair in the excited state matches exactly and, thus, 
the resonant tunneling and dissociation regions do not overlap. 
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VI. APPENDIX 



The transfer matrix propagates the plane wave solution ipi = exp(m/) from the asymp- 
totic region / < — L to the asymptotic region / > L, where the wave function is given by 
the superposition of two plane waves, ipi = aexp(m/) + 6exp(— m/). Then the tunneling 
amplitude is given by t = 1/6*. The latter relation follows from the equation 

(10) 



s 



s 



t r* 



where 5* is the scattering matrix. The relation t = 1/6* is valid independent of whether 
the wave function is scalar or vector, although the explicit form of the transfer matrix is 
different. 

We begin with the one-state model ([7]) where we assume /„, = I6m,i for simplicity: 



- - {ipi+i + ipi^i) + 2tiipi = Eipi , E{k) = -IcosK . 
It immediately follows from f ITT]) that the plane wave can be propagated as 



(11) 




Ti 



(12) 



where 



T, 



(13) 



V 



2{2ei-E)/I -1 
1 

is the one-step transfer matrix. It is worth noting that the matrix ( I13p is unitary, with the 
eigenvalues lying on the unit circle. 

The case of the two-state model is more involved. Here Eq. (fT2|) takes the form 



(1) 



(1) 



(14) 
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and the transfer matrix is given by 



Ti 




(15) 



where a = ^/2{A + + e^+i — E)/ J and h = -\/2(2e;+i — E)/ J . The dispersion relation 
E = E{k) entering equation (fT5l) is obtained by solving the Schrodinder equation in the 
absence of the scattering potential and is given in Eq.(j4l). 

Finally, we display the transfer matrix for the three-state model. It has the most compact 
form if we propagate the column vector ^P"*, ^^^^Y^ with the shifted index for the third 

component. Then the transfer matrix is given by 



Ti = A-'B , 

where the matrices A and B are written below : 

/ -I/V2 {Uo + 2ei+,~E)/J 



(16) 



A 



and 



V 



B 



-1/2 






-I/V2 

I/V2 
1/2 



\ 

(e/ + €1+2 - E)/J 
-1/2 



(17) 








A 




(18) 



y^E-Ui-ei- ei+,)/J I/V2 1/2 / 
In the paper, when discussing the tunneling probability, we focussed on the particular case 
K, = 71/2 where the bound pair has the maximal group velocity. It is interesting to compare 
results of the one-, two-, and three-state models for other values of the quasimomentum. 
This comparison is given in Fig. [7] which shows the tunneling probability as the function 
the quasimomentum k and the amplitude V of the external Gaussian potential {a = 0.65). 
The system parameters are J = 1 and Uq = Ui — 2, which implies A = 2 in the two-state 
model and / ~ 0.5 in the one-state model. (In Fig. [Tl^a) we used I = 0.7321/2, which is one 
half of the actual band width.) A narrow window of the resonant tunneling (as predicted 
by the two-state model) and partial suppression of this resonant tunneling (as predicted by 
the three-state model) are clearly seen in the figure. 
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FIG. 7: Tunneling probability as predicted by the one-state (upper left), two-state (upper right), 
and three-state (lower raw) models. The system parameters are J = 1, A = 2, and Ui = —2 in the 
lower left panel and Ui = in the lower right panel, the width of the Gaussian scattering potential 
a = 0.65. 
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